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1 Foreword 


The theta functions introduced by Jacobi [J] (see also [B . [W], [WW], [Mi]) are doubly 
(quasi)periodic analogues of the basic trigonometric functions sin(7rM) and cos(7rM). Let 
the two (quasi)periods be 1 and r G C with the condition > 0. The basic theta 
functions are 6 i{u\t), 62 {u\t), 6 ^{u\t), 6 i{u\T). The theory of theta functions is a sort of 
“elliptically deformed” trigonometry. In essence the functions sin and cos are the same 
because cosx = sin(a; + |), but everybody knows that in practice it is more convenient to 
work with the two functions rather than one. Likewise, the four theta functions can be 
obtained from any one of them by simple transformations like shifts of the argument and 
multiplying by a common factor, but it is more convenient to deal with the set of four instead 
of one. 


The “elliptic deformation” of the trigonometric functions may go in two ways depending 
on which property of the former one wants to preserve or generalize. One is a deformation 
in the class of entire functions (the north-east arrow in the diagram below). It leads to 
the quasi-periodic theta functions, which are regular functions in the whole complex plane. 
The other one is in the class of doubly periodic functions. The (inhnite) second period of 
the trigonometric functions becomes hnite (equal to r) at the price of breaking the global 
analyticity, so the elliptic functions sn, cn and dn, which are doubly periodic analogues of 
trigonometric sin and cos are meromorphic functions in the complex plane. 


{sin, cos} 


( 01 , 02 , 03 , 84} 


{sn, cn, dn} 


In fact the basic elliptic functions are constructed as ratios of the theta functions and in this 
sense the latter seem to be more fundamental. 

In practical calculations with trigonometric functions (and their hyperbolic cousins), one 
needs just a few identities for the basic functions sin and cos like the addition formula 
sin(a; + y) = sin x cos?/ J- sin?/cos x. It is not difficult to remember them all or derive any 
forgotten one from scratch using the dehnitions sinx = —?(e“—e“®^)/2, cosx = (e“-l-e““)/2. 
For the theta functions, the situation is much more involved. They are connected by a 
plethora of identities most of which are not obvious, not suitable for memorizing and can 
not be derived from scratch in any easy way. Here is what Mumford wrote in Chapter 1 of 
his book “Tata lectures on Theta I” |MI] after presenting a list of ponderous identities for 
theta functions: 


“We have listed these at such length to illustrate a key point in the theory of theta 
functions: the symmetry of the situation generates rapidly an overwhelming number 
of formulae, which do not however make a completely elementary pattern. To obtain 
a clear picture of the algebraic implications of these formulae altogether is then not 
usually easy. ” 
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All this is aggravated by the fact that there are several different systems of notation for theta 
functions in use. 

In the present paper we make an attempt to bring some order into this conglomeration of 
formulae. We show that the 3, 4 and 5-term identities of degree four (i.e. with products 
of four theta functions in each term), referred to as Weierstrass addition formulae, Jacobi 
relations, and Riemann identities, respectively, can be obtained by purely algebraic manip¬ 
ulations from six basic 3-term theta relations of degree two connecting theta functions with 
modular parameters r and 2r. Starting with the six “elementary bricks”, it is possible to 
derive 52 fundamental relations of degree four containing four independent variables. Be¬ 
sides, we give the complete list of all important particular identities which are appropriate 
specihcations of the basic bilinear and degree four ones. Recently, Koornwinder has proved 
[K] that the Weierstrass addition formulae and the Riemann identities are equivalent. We 
reproduce this result in a very simple way. 

In the future we plan to address more specihc questions related to the role of theta functions 
in the theory of integrable systems and lattice models of statistical mechanics. 
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2 Theta functions 


2.1 Theta functions with characteristics 


Fix the modular parameter r G C such that Q'r > 0 and consider the inhnite series HE]: 


da,b{u\T) = exp < 7 liT{k J- J- 2 TTi{k -|- a)(M -|- 


( 2 . 1 ) 


where i = \/—1 and a, 6 G M. The series is absolutely convergent for any u G C and dehnes 
the entire function 6 a,b{u\T). It is called the theta function with characteristics a,b. These 
functions are connected by the relations 

OaAu + ra' + b'\T) = 

0a+l,b{u\T) = 9a,b{u\T), 

6a,b+i{u\r) = e^™6'a,fe(M|r). 
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(2.2a) 

(2.2b) 

(2.2c) 




In particular, it follows from here that the functions 6*a ^ are quasiperiodic with (quasi)periods 
1 and r: 


9a,b{u + l|r) = e^'^"°‘9a,b{u\T), (2.3a) 

9aAu + r|r) = ( 2 . 3 b) 


and the shifts by half-periods are given by 

9a,h{'^ + \\t) = 0„,b+i(M|r) = (M|r), (2.4a) 

Oa,b{u + i|r) = (2.4b) 

It follows from the definition that 9a,b{—u) = 9_a-b{u). Hence, according to (I2.2bp . fl2.2cp . 
the functions 9a,b{u) have definite evenness properties only for integer or half-integer char¬ 
acteristics. By virtue of fl2.2bp . fl2.2cp . it is sufficient to consider the theta functions with 
characteristics 0 < a, 6 < 1. 


2.2 Basic theta functions 

The standard theta functions with half-integer characteristics mm are defined as follows: 
9i{u\t) = 

92{u\t) = 

93{u\t) = 

94{u\t) = 

where 

q := e "*", | g | < 1 . ( 2 . 6 ) 

19 19 

In the limit r —)■ ioo they are: 9i{u\t) = 2qi sinvrM -|- 0{qi), 92 {u\t) = 2qi costtm -|- 0{q'^), 
9i,{u\t) = 1 -h 0{q), 94 {u\t) = 1 -F 0{q). 

In what follows we often write 9r{u\T) := 9r{u), r = 1, 2, 3,4 if this does not cause confusion. 
From fl2.5p it is clear that the function 9i is odd, 9i[—u) = — 6 'i(m); the other three are even, 
9s{-u) = 9s{u), s = 2,3,4. 

The values 6*((0), 6 * 2 ( 0 ), 6 * 3 ( 0 ), 6 * 4 ( 0 ) are called theta constants. 

2.3 Shifts by periods and half-periods 

Here we list the essential transformation properties for the theta functions fl2.5p which follow 
from fl 2 . 2 p . 


k ■Ki{2k+l)u 


-9i^i{u\t) = -i^{-iyq^'" e 

kei, 

e,yu\T)= 


fcez 

6*0,0 (m k) = 


27riku 


(2.5) 




k k^ 27Tiku 

q 6 , 
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Shifts by periods: 


9i{u + 1) = -di{u), 
02 (m + 1 ) = -02 (m), 
03(m + 1) = 03(m), 

04(m + 1 ) = 04(m). 


0i(M + r) = -e-"(2«+")0i(M), 

02(M + r) = 

04(M + r) = -e-"(2“+")04(M). 


0i(M + r+l) 

02 (M + r+l) = -e-^'( 2 «+") 02 (M), 
03(M + r+l)=e-"(2“+")03(M), 
04 (M + r+l) = 


Shifts by half-periods: 


01 (m + |) = 02 (m), 
d2{u+ \) = -01 (m), 
03(m + i) = 6'4 (m), 
04(m+ \) = 03(m). 


0i(M + f) =ie-"'(“+"/^)04(M), 
02 (M + i) = e-^'(“+"/^)03(«), 
03(m + i) = e-"'(“+"/4)02(M), 

04(m + i) =ie-”(“+-/4)04(M). 


01 («+^) 
02 («+^) 
03(«+^) 


_ p—-Ki(u+Tl 


e ''^^“"/^^ 03 (m), 

-"(“+"/4)04(m), 

ie-"'(“+"/4)0i(M), 


= —*e 


_ ^-ni{u+T 




(2.7) 


( 2 . 8 ) 


2.4 Zeros of theta functions 


These relations imply that the (first order) zeros of the theta functions are as follows: 


where n, m G Z. 


01 (n) = 0 
02 (n) = 0 
03 (n) = 0 
04 (n) = 0 


M = n -|- mr, 
u = n + \ + mr, 
u = n + \ + {m + \)t, 
u = n + {m+ |)r, 


(2.9) 


Indeed, in accordance with (12.7p . (|2.8p . the functions 0i(n-|-n-|-mr), 02(rt + n-|-^-|-mr), 9^{u+n+ 
^-|-(m-|-^)r), 01 (M-|-n-|-(m-|-^)r) are proportional to the odd function 0 i(m). Hence the corresponding 
zeros are as in p2.9p . To complete the proof, it is sufficient to show that the function 0 i(m) has 
precisely one simple zero in the parallelogram H with the vertices — ^ ± ^^ ± ^. The standard 
argument is to compute the contour integral ^ fgjj ^ (^) = f which means that the zero is 

simple (see [WW] . [Ml] for details). 
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2.5 Theta functions as infinite products 


One has the following inhnite product representations: 

CXD 

ei{u\T) =2g3sin7ruJJ(l-g2«)(l-g2V™)(l-g2’^e-2™), (2.10a) 

n=l 

oo 

02(M|r) =2g3cos7rM JJ(l-g2")(l + g2V™)(l + g2"e-2™), (2.10b) 

n=\ 

oo 

03 (M|r) = JJ(1 - g2-)(l + ^ ^2n-l ^ ^2.10c) 

n=l 

oo 

e^{u\T) = JJ(1 - g2-)(l _ g2n-lg2™^^^ _ ^2n-lg-2™^_ ^2.10d) 

n=l 

To prove (|2.10dp . we note that the product p{u\t) := n^i(l “ q^"'“^e^™)(l — g2n - ig - 2 ™^ 
the same zeros as 04 (u|r) and the ratio 04 (ri|r)/p(ri|r) is a doubly periodic function with periods 
1 and T. Hence the ratio is constant and one has 04 (ri|r) = Ap{u\t). To find the constant H, put 
u = 0 thus getting A = 04 (O|r)/p(O|T). Finally, in accordance with the Gauss formula |Anl p. 23, 
eq. (2.2.12)], 

1 Th 

0,m ==n (2-11) 

keZ n=l ^ 

Rewriting n^i(f “ ~ q^^)i oil® gets A = n^i(l “ 9^”) and formula 

(j2.10dp is proved. Equations (I2.10ap - (|2.10cp are obtained from (I2.10cp by appropriate shifts of u 
in accordance with (12.81) . 

As a corollary of fl2.10p one has inhnite product representations for the theta constants: 


e[{0) = 27rq-^llil-q^^f, (2.12a) 

n=l 

oo 

02(0) = 2g2 JJ(1 - g2-)(l + g2n)2^ (2.12b) 

n=l 

oo 

^ 3 ( 0 ) = Hil - g'")(l + (2.12c) 

n=l 

00 

04(0) = JJ(1 - g2^)(l - g2«-i)2. (2.12d) 

n=l 


Since ]^(1 + 0'^”)(1 + g^"'~^)(l — g^"'“^) = 1, this implies the famous identity for the theta 

n>l 

constants HI p. 517]: 


0'l(O) = 7r02(O)03(O)04(O). 


(2.13) 
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2.6 Modular transformations 


The transformation r ^ r + 1: 


0i{u\t + 1) = e 4 01 (u r). 

(2.14a) 

02 (m t + 1) = e^02(M r). 

(2.14b) 

03(M|r + 1) = 04(^1^), 

(2.14c) 

0/^{u\t + 1) = 03 (ur). 

(2.14d) 

Since r ^ r + 1 implies q —> —q. equations fl2.14p follow from 02.51) or 02.1()p. 

The transformation r — I/'t: 


01 {u/t\ — 1/r) = 

(2.15a) 

02 {u/t\ — 1/r) = -s/— ire™ ^'”04 (m r). 

(2.15b) 

03 (u/r - 1/r) = \/-ire™‘'/^03(M r). 

(2.15c) 

04 {u/t\ — 1/r) = V— iT e™*^'”02(M r). 

(2.15d) 


The branch of the square root here is such that 'R\/—iT > 0. 

The proof is well known. Since the ratio {u/t\ — 1/r) / 03 (M|r) := C is an entire doubly 

periodic function of u, it is a constant (which may depend only on r). Shifting u by one 

obtains three more relations of the same kind with the same constant C. Then the substitution 
of these formulas to ()2.13p . yields = —ir. The sign of the square root is determined by the 
argument that if r G then both 03 ( 0 | — 1/r) and 03 (O|r) are real and positive. 


2.7 Other notation for the theta functions 


The notations for theta functions used in the literature are of a great variety. This can be 
a source of confusion. Here we briefly comment on the main systems of notation other that 
the one adopted in this paper. In the theory of elliptic integrals, the theta functions 


Qr{u\T) = Or 



t), K = ^01(0\t) 


( 2 . 16 ) 


introduced by Riemann are commonly encountered. The number K is the full elliptic integral 
(of the hrst kind). In [A] and some other books our Or is denoted as dr while 0^ dehned 
in (12.161) is just Or- The antiquated Jacobi notation (still preferred by some authors) are if, 
ifi, © 1 , 0 for 01, 02, 03 , 04 respectively. The “multiplicative notation” 0r{z\q) for 0j,(w|r), 
where q = , z = is widely used in the modern literature on elliptic hypergeometric 

series and related problems. 


Lastly, let us mention a few of the minor differences in notation encountered in the literature. 
In [W] , [HCj the functions Qa,b{.u) have been considered which are related with Oa^iu) by 


eUu) = e--^0.dul QZ(u) = e 


TTzab 

2 9 a h\ 
2 ’ 2 


u] 


(2.17) 
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The set of our theta functions fl2.5l) is related with the corresponding functions in [WW] as 
6r{u\T) = r = 1,2, 3,4. Following the original notation [J], in and in some 

other books the notation 6q is used instead of 6 * 4 . 


3 Four types of identities between theta functions 

3.1 Preliminaries 

The number of identities satished by the theta functions is enormous. It is still fairly big 
if we consider identities involving up to four independent variables. They can be split into 
four types: 

B. Three-term bilinear identities involving two independent variables. They relate 
products of two theta functions with modular parameter r to linear combinations 
(actually, sums or differences) of similar products of theta functions with modular 
parameter 2r. 

W. Three-term identities of degree 4 (the Weierstrass addition formulae). 

J. Four-term identities of degree 4 (the Jacobi formulae). 

R. Five-term identities of degree 4 (the Riemann identities). 

The identities of types W, J, R include theta functions with the same modular parameter r 
and contain four independent variables. The identities of type B are the most fundamental 
ones: all the others are algebraic consequences of these together with the evenness properties 
of the theta functions 9r{—u) = {—l)^^'^9r{u), r = 1,2, 3,4. Namely, we shall show how to 
derive W from B etc., according to the scheme B — )■ W — )■ J —)■ R. It also turns out that 
the Jacobi and Riemann identities are equivalent in a very simple way. At the end of this 
section, we prove the arrow W 4 — J which implies equivalence of the Weierstrass and Jacobi 
identities. 


3.2 Three-term bilinear identities connecting theta functions with 
r and 2r 


B.I. There are six basis bilinear identities: 


9i 

{u 

t)9i( 

V 


= 93 

[u 

+ 

V 

2t)92 

>- 

- V 

9i 

{u 

t)92 

[v 


= 93 

[u 

+ 

V 

2r)04 

[u- 

- V 

92 

{u 

t)92 

[v 


= ^2 

[u 

+ 

V 

2r)03 

[u- 

- V 

93 

{u 

r)93 

V 


= 93 

[u 

+ 

V 

2r)03 

[u- 

- V 

O 3 

{u 

t)94{ 

V 


= 04 

[u 

+ 

V 

2r)04 

[u- 

- V 

9, 

{u 

t)94{ 

V 


= ^3 

[u 

+ 

V 

2r)03 

>- 

- V 


2t) — 92{u + v\2t)93{u — v\2t), (3.1a) 

2t) + 9i{u + v\2t)9i{u — v\2t), (3.1b) 

2t) + 93{u + v\2t)92{u — v\2t), (3.1c) 

2t) + 92{u + v\2t)92{u — v\2t), (3.Id) 

2r) — 6 'i(m J- n|2r)6'i(M — n|2r), (3.1e) 

2t) — 92{u + v\2t)92{u — v\2t). (3.If) 










(See Ig , [D], |M2] for the general case of multi-dimensional theta functions.) 
B.II. A system equivalent to fl3.ip : 

29i{u + v\2t)9i{u — v\2t) = 9i{u\T)93{v\T) 

01 (u + v\2t)9a(u — v\2t) = 01 (M|r) 09 (n|r) 


-,-w ,• / 6 ' 3 (M|r) 04 (n|r), 

'i(m -|- n| 2 r) 04 (-u — n|2r) = 0 i(M|r) 02 (n|r) -|- 02 (n|r) 


20 .,- , 

202 (m + n| 2 r) 02 (M - n|2r) = 93{u\t)93{v\t) — 04 ( 
202 (m + n| 2 r) 03 (M - n|2r) = 02 (M|r) 02 (- 
29^{u + n| 2 r) 03 (M — n|2r) = 9, 


9'a(u\t)9'a( 


>|r) 04 (n|r), 
n|r) - 0i(M|r)0i(n|r), 


(3.2a) 

(3.2b) 

(3.2c) 

(3.2d) 

(3.2e) 

(3.2f) 


Remark 3.1 Starting with any identity in \3.1\] . one ean derive all the other ones by ap¬ 
propriate shifts of the variables u, v. 

The proof is standard. Let us prove, for example, (I3.1bp . Consider the function 

F{v) := 6i(u -|- n|2r)04(n — n|2r) -|- 04 (m -|- n|2r)0i(u — n|2r). 

By virtue of (j2.7p and (j2.8|) . F{v -|- 1) = —F{v), F{v + t) = F{v) and F{^) = 0. 

Hence zeros of F{v) are Vn,m = n -\- ^^ + n^'r, n,m ^ X and the ratio F{v)/{ 9 i{u\t) 62 {v\t)) is an 
entire function doubly periodic in v with periods l,r. Therefore, this ratio does not depend on v: 
F{v)/ 9 i{u\t) 62 {v\t) = C{u). Setting v = u, one has: 

^ 0 i( 2 n| 2 r) 04 (O| 2 r) 

0 i(n|T) 92{u\t) 

By virtne of (I2.10al] . (|2.10bp and (|2.1ip . 9 i{u\t)92{u\t) = 0 i( 2 n| 2 r) 04 (O| 2 r) and thus C{u) = 1. 


3.3 Three-term Weierstrass addition identities 


There are twelve addition formulae (see below). We start with the identity 


01 (u -I- a;)0i(M — x)9i{v + y)9i{v — y) — 0i(m + y)9i{u — y)9i{v + x)9i{v — x) 

= 9i{u-\- v)9i(u — v)9i{x y)9i{x — y) 


(3.3) 


which was originally discovered and proved by Weierstrass [wa p. 155]. All the identities 
listed below in this section can be derived from it by appropriate shifts of the variables in 
accordance with relations 02.81) . Our approach is different. We show that all the identities of 
Weierstrass’ type are simple algebraic consequences of the bilinear system 03. ip together with 
the evenness conditions 9r{—u) = {—lY''’^9r{u), r = 1,2, 3,4. This argument is independent 
of dUp. 

To prove (|3.3I) . one should rewrite (|3.1al) as 


6 i{u x \ t ) 6 i{u — x \ t ) = 03(2M|2r)02(2a:|2r) — 02(2u|2r)03(2x|2r). 
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Multiply this by the similar expression for 6 i{v + y|r)0i(u — y|r) and subtract the same with the 
change x o y. Using (I3.1al) once again, we arrive at (13.3p . All the equations below in this section 
can be obtained from system (|3.ip in a similar way. 

W.I. Symmetric system: 


Oiiu + x)6i{u — x)6r{y + y)0r{v — y) — 6i{y + x)6i{y — x)6r{u + y)6r{u — y) 
= 6i{u + v)6i{u — v)9r{x + y)0r{x — y), 


(3.4) 


r = l,2,3.4. 

W.II. Complimentary system: 


92{u + x)92{u - x)93{v + y)93{v - y) - 92{v + x)92{v - x)9‘i{u + y)9z{u - y) 
= —9i{u + v)9i{u — v)9i{x + y)9i{x — y), 


(3.5a) 


92{u + x)92{u - x)9i{v + y)9i{v - y) - 92{v + x)92{v - x)9i{u + y)9i{u - y) 
= -9i{u + v)9i{u - v)9^{x + y)9^{x - y), 


(3.5b) 


92,{u + x)92,{u - x)9i{v + y)9i{v - y) - 9^{v + x)9^{v - x)9i{u + y)9i{u - y) 
= -9i{u + v)9i{u - v)92(x + y)92(x - y). 


(3.5c) 


W.III. Asymmetric system: 

9r{u + x)9r{u — x)9r{v + y)9r{v — y) — 9r{u + y)9r{u — y)9r{v + x)9r{v — x) 
= {-iy~^9i{u + v)9i{u - v)9i{x + y)9i{x - y), 

r = l,2,3,4. 

W.IV. Complimentary identity: 

9yu + x)9yu - x)9yv + y)9yv - y) - 9yv + x)9yv - x)9yu + y)9yu - y) 
= 02 (m + v)92{u - v)92{x + y)92{x - y). 


(3.6) 


(3.7) 


W.V. Mixed identity: 

01 (m + x ) 92 {u - x)9yv + y)04(f - y) - 9i{u - y)02(M + y)9'i{v - x)9yv + x) 
= 01 (x + y) 02 (a; - y) 03 (M + v)9yu - v). 


(3.8) 


Remark 3.2 Sometimes the Weierstrass addition formula LS. ,9\) is referred to as Fay iden¬ 
tity. In fact, it is a generalization of Jacobi’s results (see Section [7^ below). 
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3.4 Four-term Jacobi identities 


In Sect. 13.31 we have presented twelve three-term identities of degree four depending on four 
variables u, v, x, y. Here we introduce another set of variables IH, X, F, Z and their “dual” 
counterparts [WW] : 

W' = 1{-W + X + Y + Z), 

X'=l{W-X + Y + Z), 

? (3.9) 

Y' = ^{W + X-Y + Z), 

Z' = 1{W + X + Y-Z). 

One can easily verify that W,X,Y,Z are expressed via the “dual” variables W', X', Y', Z' 
by the same formulae, i.e., the “prime procedure” applied to fl3.9p yields (IF')' = W = 
|(—hh' -|- X' -|- F' -|- Z') etc. We employ the short-hand notation 

\pqrs\ := 0p(lF)0,(X)0,(F)0,(Z), \pqrs\ := 0p(lF')0,(X')0,(F')0,(Z') 

which is widely used in [WW] . If all the indices of the theta functions coincide, this is 
further abbreviated to 

[r] ;= 0,(1F)0,(X)0,(F)0,(Z), [r]' := 0,(1F')0,(X')0,(F')0,(Z'). 

Below we list all four-term basic identities of degree four which were essentially obtained by 
Jacobi [J] p. 507]. Here we present these in a more symmetric and comprehensive form. 

The simplest (and most important) ones are: 



[1] + [2] = [1]' + [2]', 

(3.10a) 

J.I 

[1] - [2] = [4]' - [3]', 

(3.10b) 


[3] + [4] = [3]' + [4]', 

(3.10c) 


[3] - [4] = [2]' - [1]'. 

(3.10d) 


The system (I3.10p is a direct algebraic corollary of appropriate addition formulae given in Section 
[T3l To see this, we relate the variables u, v, x, y with the variables of the present section as 
follows: 


W = u + X, 
X = u — X , 

< 

Y = v + y, 

_ Z = v-y. 


' W' = v-x, 
X' = V + X , 
Y' = u-y, 

^ Z' = u + y. 


(3.11) 


Further, the products of theta functions containing “inappropriate” combinations tt± u, x±y can 
be excluded from the addition formulae. Then identities (I3.ina|] . (jJ.lOcp emerge as particular cases 
of (|3.6p . Changing u x in ()3.6p (with r = 3) and ()3.7p . one obtains (|3.10bh . Finally, (|3.10dp is a 
“dual” version of (|3.10bp . 
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Remark 3.3 Equations hS-Kh) are a pari of the system of twelve identities written in [WWl 
pp. 468, 488]. It is easy to see that all additional relations are appropriate linear combina- 


the basic ones. Hd.l(J\). 

For completeness, 

we give here the full list: 

[1] + 

[2] = 

= [1]' + [2]', 

[1] 

+ [3] 

= [2]' + 

[4]', 

[1] 

+ [4] = 

[1]' + [4]', 

[2] + 

[3] = 

= [2]' + [3]', 

[2] 

+ [4] 

= [1]' + 

[3]', 

[3] 

+ [4] = 

[3]' + [4]', 

[1]- 

[2] = 

= [4]' - [3]', 

[1] 

-[3] 

= [1]' - 

[3]', 

[1] 

-[4] = 

[2]' - [3]', 

[2]- 

[3] = 

= [1]' - [4]', 

[2] 

-[4] 

= [2]' - 

[4]', 

[3] 

-[4] = 

[2]' - [1]'. 


(3.12) 


Now we list symmetric self-dual identities for products of type [rrss] which can also be 
derived from the addition formulae by algebraic manipulations: 


[1122] + [2211] = [1122]' + [2211]', (3.13a) 

[1133] + [3311] = [1133]' + [3311]', (3.13b) 

J.II [1144] + [4411] = [1144]' + [4411]', (3.13c) 

[2233] + [3322] = [2233]' + [3322]', (3.13d) 

[2244] + [4422] = [2244]' + [4422]', (3.13e) 

[3344] + [4433] = [3344]' + [4433]'. (3.13f) 

Farther, there are simple complimentary relations: 

[1122] - [2211] = [3344]' - [4433]', (3.14a) 

[1133] - [3311] = [2244]' - [4422]', (3.14b) 

J.III [1144] - [4411] = [2233]' - [3322]', (3.14c) 

[2233] - [3322] = [1144]' - [4411]', (3.14d) 

[2244] - [4422] = [1133]' - [3311]', (3.14e) 

[3344] - [4433] = [1122]' - [2211]'. (3.14f) 


The snbsystems (I3.14al) - (l3.14cj) and fl3.14dl) - fl3.14fll are dual to each other. Thus, the systems 
(I3.13p . (13.141) can be represented in very compact form: 

[rrss] -f [ssrr] = [rrss]' -f [ssrr]', 

[rrss] — [ssrr] = [rrss]' — [ssrr]', 

where r, s G (1, 2, 3,4), r < s and s, r G (1, 2, 3,4)\(r, s), r < s. 

Finally, there are fonr “fnlly mixed” identities: 

[1234] + [2143] = [3412]' + [4321]', 

J.IV [1234] - [2143] = [2143]' - [1234]', 

[3412] + [4321] = [1234]' + [2143]', 

[3412] - [4321] = [4321]' - [3412]'. 


(3.15) 


(3.16a) 

(3.16b) 

(3.16c) 

(3.16d) 
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Identities (|3.13al) - (l3.13cl) follow from (|3.4p . Indeed, one can write (13.41) as 


[llrr] — [llrr]' = 6i{u + v)9i{u — v)6r{x + y)0r{x — y), r = 1, 2,3,4. 

Changing here a: -H- y, one gets [rrll]'—[rrll] = 9i{u+v)9i{u—v)6r{x+y)9r{x—y) = [llrr] —[llrr]'. 
Similarly, (|3.13dp - (l3.13fp follow from (|3.5ap - (l3.5cp . respectively. To prove (|3.14al) . we write (I3.13al) 
in terms of the variables u, v, x, y and exchange u x. Then 

[2211] — [1122] = —9i{u + v)9i{u — v)92{x + y)92(x — y) + 92{u + v)92{u — v)9i{x + y)9i{x — y). 

Now (I3.14ap holds by virtue of ()3.5cp . Identities (|3.14bp - (|3.14fll can be proved in a similar way. 
Finally, it is easy to see that (j3.16p follows from (13.8p . Indeed, subtracting (|3.8p from the same 
identity with the exchange x y yields (I3.16ap . All other identities in (13.161) are proved in a similar 
way. 


Remark 3.4 Identities ^3.13\} . ^3.14\ ) differ slightly from those written by Jacobi. For 
example, in [H p. 507] one can find the relations [1122] + [4433] = [2211]' + [3344]', 
[1122] — [4433] = [1122]' — [4433]' which are appropriate linear combinations of 1,9. 13a\} . 
^3.13j\ ), Hid.lfal ), i3.14j\ ). We should also stress that these Jacobi identities are direct corol¬ 
laries of i3.4\) at r = 2 and 113. Sd) . This is in complete agreement with the derivation of 
h3.13t) . \3.14\) from the Weierstrass addition formulae. 


3.5 Five term Riemann identities 

The Riemann identities (the term is due to Mumford [Mil page 20]) are simple corollaries 
of the Jacobi relations fl3.10p . fl3.13p - fl3.16p . They each express a “primed” quantity as a 
linear combination of some appropriate four “unprimed” ones. 

Hence from (13.101) we have the four simplest Riemann identities: 



2|1|' = |1| + |2| - |3| + |4|, 

(3.17a) 


2|2]' = |1] + |2] + |3] - |4], 

(3.17b) 

R.I 

2|3]' = -11] + [2] + |3| + |4|, 

(3.17c) 


2|41' = |1] - [2] + |3] + |4], 

(3.17d) 


Let us emphasize that fl3.17p is equivalent to (13.101) . 

Next we list all possible (twelve) identities that are obtained from fl3.13p . fl3.14|) : 

2[1122]' = [1122] + [2211] + [3344] - [4433], (3.18a) 

2[1133]' = [1133] + [3311] + [2244] - [4422], (3.18b) 

2[1144]' = [1144] + [4411] + [2233] - [3322]. (3.18c) 
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2[2211]' = [2211] + [1122] + [4433] - [3344], (3.19a) 

2[2233]' = [2233] + [3322] + [1144] - [4411], (3.19b) 

R.II 2[2244]' = [2244] + [4422] + [1133] - [3311]. (3.19c) 

2[3311]' = [3311] + [1133] + [4422] - [2244], (3.20a) 

2[3322]' = [3322] + [2233] + [4411] - [1144], (3.20b) 

2[3344]' = [3344] + [4433] + [1122] - [2211]. (3.20c) 

2[4411]' = [4411] + [1144] + [3322] - [2233], (3.21a) 

2[4422]' = [4422] + [2244] + [3311] - [1133], (3.21b) 

2[4433]' = [4433] + [3344] + [2211] - [1122]. (3.21c) 


Clearly, the system (I3.18p - (l3.2ip is equivalent to (I3.13p . (I3.14p . Finally, identities (I3.16p are 
equivalent to 

2[1234]' = -[1234] + [2143] + [3421] + [4312], (3.22a) 

2[2143]' = -[2143] + [1234] + [3412] + [4321], (3.22b) 

R.III 2[3412]'=-[3412] + [4321] + [1234] + [2143], (3.22c) 

2[4321]' = -[4321] + [3412] + [1234] + [2143]. (3.22d) 

Remark 3.5 The identities presented here essentially coincide with the ones given by Mum- 
ford |M1[ p. 20]. See also [WW] . 

3.6 Equivalence of addition formulae and Jacobi identities 

In section 1331 we have obtained the Jacobi identities from the addition formulae. In its turn, 
one can show that the system (I3.10p . (I3.13p . (I3.14p implies the addition formulae (I3.5p - (l3.8p . 
The proof is similar to the one given by Koornwinder for the Riemann identities [K] . 

In accordance with (13.lip , the relation (I3.10al) acquires the form 

9i{u + x)9i{u - x)9i{v + y)9i{v - y) + 92{u + x)92{u - x)92{v + y)92{v - y) 

(3.23) 

= 9i{v + x)9i{v - x)9i{u + y)9i{u - y) + 92{v + x)92{v - x)92{u + y)92{u - y). 

Changing here u ^ x and n -H- x, one obtains two additional relations: 

-9i{u + x)9i{u - x)9i{v + y)9i{v - y) + 92{u + x)92{u - x)92{v + y)92{v - y) 

(3.24) 

= -9i{u + v)9i{u - v)9i{x + y)9i{x - y) + 92{u + v)92{u - v)92{x + y)92{x - y). 
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(3.25) 


9i{u + v)9i{u - v)9i{x + y)9i{x - y) + 92{u + v)92{u - v)92{x + y)92{x - y) 

= -9i{v + x)9i{v - x)9i{u + y)9i{u - y) + 92{v + x)92(v - x)92{u + y)92{u - y). 

Introduce the notation: 

Aj := 9j{u + x)9j{u — x)9j{v + y)9j{v - y), 

~ + x)9j{v — x), (3.26) 

Cj := 9j{u + v)9j{u — v)9j{x + y)9j{x — y). 

Then relations (I3.23p - (|3.25p acquire the form Ai — Bi = B 2 — A 2 , Ai — Ci = A 2 — C 2 , Bi + Ci = 
B 2 — C 2 which is a system of linear equations for the unknowns A 2 , B 2 ,C 2 - The system is degenerate 
with compatibility condition Ai — Bi = Ci. In terms of the theta functions, this condition is nothing 
but equation (13.6p with r = 1. Equivalently, one can treat the above equations for Aj,Bj,Cj as a 
linear system for the unknowns Ai,Bi,Ci. Then, for example, Ci = B 2 — A 2 which is (13.6p with 
r = 2. The other addition formulae can be obtained in a similar way. 


4 Particular identities 


4.1 Consequences of the bilinear identities 

One can obtain twelve particnlar identities from the general system (13.2p putting n = 0 or 
V = ±u (actually, the restriction v = —u can be applied only for identity fl3.2dp which leads 
to (iroil below): 

2el{u\2T) = 64 {u\t) 63 { 0 \t) - 6 ' 3 (M|r) 6 ' 4 ( 0 |r), (4.1a) 

2el{u\2T) = 6'3(M|r)6'3(0|r) - 6'4(M|r)6'4(0|r), (4.1b) 

291{u\2t) = 6^3 (m I r) 6 ^ 3 ( 0 1 r) + 6 ' 4 (M|r) 6 ' 4 ( 0 |r), (4.1c) 

26'4(M|2r) = 6'3(M|r)6'4(0|r) + 6'4(M|r)6'3(0|r), (4.Id) 

26*i(M|2r)6'4(M|2r) = 6 'i(M|r) 6 * 2 ( 0 |r), (4.2a) 

26'2(M|2r)6'3(M|2r) = 6'2(u|r)6'2(0|r), (4.2b) 

26*2 (2m 12 r) 6*2(0 12 r) = (^3{u\t) — 6*4(M|r), (4.3a) 

26*2 (2m 1 2 r)6*3(0 1 2 r) = ^2('*^k) ~ (4.3b) 

203 ( 2 M| 2 r) 02 (O| 2 r) = el{u\T) + 0?(M|r), (4.3c) 

203(2M|2r)03(O|2r) = 91{u\t) + 6 * 4 (M|r), (4.3d) 

6*1 (2m 1 2r) 6 * 4(01 2r) = 6 *i(M|r) 6 * 2 (M|r), (4.4a) 

6*4 (2m 1 2r) 6 * 4(0 1 2r) = 6 * 3 (M|r) 6 * 4 (M|r). (4.4b) 
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In |WW[ Section 21.52] there are two particular equations relating theta functions with 
modular parameters r and 2r which are called the transformations of Landen’s type: 


94 { 2 u \ 2 t ) 6*3(M|r)6*4(M|r) 


04 (O| 2 r) 03 (O|r) 04 (O|r)’ 

9i{2u\2t) 6*i(M|r)6*2(M|r) 


04 (O| 2 r) 03 (O|r) 04 (O|r)- 

The hrst identity is derived from fl4.4b[) and from the relation 

0l(O|2r) = 03(O|r)04(O|r) 


(4.5a) 


(4.5b) 


(4.6) 


which is also a corollary of fl4.4bp . The identity fl4.5bp is a ratio of fl4.4ap and fl4.6p . 
Note also that fl4.5bp can be obtained from fl4.5ap by the shift m —)■ k + |. 


4.2 Particular addition formulae 

One can obtain important particular cases of the Weierstrass addition formulae which include 
two variables. Here we present the complete list of eighteen identities which easily follow 
from fl3.4p - fl3.8p : 


9i{u + v)9i{u — v)9l{0) = 9l{u)9l{v) — 9l{u)9l{v) 

= 9l{u)9l{v)-9l{u)9l{v), 


(4.7a) 


9i{u + v)9i{u — v)9l{0) = 9l{u)9l{v) — 9l{u)9l{v) 

= 9l{u)9l{v)-9l{u)9l{v), 


(4.7b) 


9i{u + v)9i{u - v)9l{0) = 9l{u)9l{v) - 9l{u)9l{v) 

= 9l{u)9l{v) - 9l{u)9l{vf 


(4.7c) 


92{u + v)92{u - v)9l{0) = 9l{u)9l{v) - 9l{u)9l{v) 

= 9l{u)9l{v)-9l{u)9l{v), 


(4.8a) 


92iu + v)92iu - v)9l{0) = 9l{u)9l{v) - 9l{u)9l{v) 

= 9l{u)9l{v) - 9l{u)9l{vf 


(4.8b) 


92iu + v)92iu - v)9l{0) = 9l{u)92{v) - 9j{u)9l{v) 

= 9i{u)9!{v)-9i{u)9!{v), 


(4.8c) 


16 

























9-i{u + v)9^{u - v)9l{0) 


= 9l{u)9l{v) + 9l{u)9l{v) 

= 9l{u)9l{v) + 9liu)9l{v), 


(4.9a) 


93{u + v)93{u - v)6'3(0) 


9l{u)9l{v) + 9liu)9l{v) 

9l{u)9l{v) + 9liu)9l{v), 


(4.9b) 


93{u + v)93{u - v)9l{0) 



(4.9c) 


9i(u + v)94{u — v)9l{0) 


9l{u)9l{v) + 9l{u)9l{v) 

9liu)9liv) + 9l{u)9liv), 


(4.10a) 


6*4(m + v)94{u — w)6*3(0) 


= 9l{u)9i{v) + 9l{u)9l{v) 

= 9l{u)9l{v) + 9l{u)9l{v), 


(4.10b) 


6*4 (m + v)94{u — v)9‘l{9i) 



(4.10c) 


9i{u + v ) 92 {u - v)6*3(0)6'4(0) = 9 i { u ) 92 { u ) 93 { v ) 94 { v ) + 93 { u ) 94 { u ) 9i { v ) 92 { v ), (4.11a) 

9i{u + v ) 93 {u - v)6'2(0)6'4(0) = 9 i { u ) 93 { u ) 92 { v ) 94 { v ) + 6'2(M)6'4(M)6'i(u)6'3('y), (4.11b) 

9i{u + v ) 94 {u - v)02(O)03(O) = 6*i(M)6*4(M)6*2(f)6'3(t>) + 92 { u ) 93 { u ) 9 i { v ) 94 { v ), (4.11c) 

6*2(m + v ) 93 {u - f)6*2(0)6*3(0) = 6*2(M)6*3(M)6*2(f)6*3(i;) - 6'i(M)6'4(M)6*i(t>)6*4(t>), (4.11d) 

6*2(m + v ) 94 {u - f)6*2(0)6*4(0) = 6*2(M)6*4(M)6*2(f)6*4(i;) - 6*i(M)6*3(M)6*i(t>)6*3(w), (4.lie) 

93 {u + v ) 94 {u - v)6'3(0)6'4(0) = 93 { u ) 94 { u ) 93 { v ) 94 { v ) - 9i { u ) 92 { u ) 9i { v ) 92 { v ). (4.11f) 


Remark 4.1 The complete list of identities fjf. 7\) - ^.ll\ ) was originally obtained by Jacobi 
[H p. 510] as a particular specification of identities 113. 10\) . ^3.13\} - l[3.1&) . see also |Wl pp. 
76-78], |WW[ 487-488]. Mumford |Mll p. 22] has obtained a part of relations (^T?)-(^T7I]) 
as specific cases of the Riemann identities \3.11\)-[3. 22i). 
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As a byproduct of fl4.7l) - fl4.10p . one gets some extra identities: 

el{u)el{v) - el{u)el{v) = el{u)9l{v) - el{u)el{v), (4.i2a) 

el{u)el{v) - 9l{u)9l{v) = 9l{u)9l{v) - 9l{u)9l{v), (4.12b) 

9l{u)9l{v) - 9l{u)9l{v) = 9l{u)9l{v) - 9l{u)9l{v), (4.12c) 

9l{u)9l{v) - 9l{u)9l{v) = 9l{u)9l{v) - 9l{u)9l{v). (4.12d) 

In particular, the following identity holds: 

9{{u) + 9l{u) = 9^{u) + 9^{u). (4.13) 

Certainly, the last relation is a corollary of fl4.3p . Thus, in addition to fl2.13l) . one gets 
another famous identity for theta constants: 

03^(0) = 02^(0)+ 01(0). (4.14) 


4.3 Duplication formulae 

The duplication formulae relate the functions 9a{2u\T), a = 1,2, 3,4 with appropriate com¬ 
binations of the functions 9b{u\T). All these identities emerge as further degenerations of 


addition formulae fl4.7ll - fl4.lip . Here is the complete list: 

0i(2n)02(O)03(O)04(O) =2 9^{u)92{u)9s{u)9,{u), (4.15) 

02 ( 2 n) 02 (O) 0 ^(O) = 9Uu)9l{u) - 9^,(u)9l(u), (4.16a) 

02(2n)02(O)0l(O) = 9l{u)9l{u) - 9j{u)9l{u), (4.16b) 

92{2u)9l{0) = 0l(n) — 9^{u), (4.16c) 

92{2u)9l{0) = 9t{u) - 9l{u). (4.16d) 

03 ( 2 n) 03 (O) 02 '(O) = 9l{u)9l{u) + 0?(n)0^(n), (4.17a) 

9s{2u)9s{0)9l{0) = 9l{u)9l{u) - 9l{u)9l{u), (4.17b) 

03 ( 2 n) 0 |(O)= 9t{u) + 9l{u), (4.17c) 

9s{2u)9l{0) = 0l(n) + 9t{u). (4.17d) 

04 ( 2 m) 04 (O) 02 (O) = 02 ('^)^ 4 ('*^) + ^l('*^)^ 3 ('^)) (4.18a) 

04 ( 2 ^) 04 ( 0 ) 03 ( 0 ) = 9l{u)9l{u) + 9l{u)9l{u), (4.18b) 

04 ( 2 n) 0 l(O) = 9l{u) - 0l(n), (4.18c) 

04(2n)0l(O) = 0l(n) - 01 (m). (4.18d) 
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One can unify the sets of equations { fl4.16ap . fl4.17bD . fl4.18al) }: { fl4.16bp . fl4.17ap . fl4.18bp }: 
{ fl4.16cp . fl4.17cp . fl4.18cp } and { fl4.16dp . fl4.17dp . fl4.18dp } by writing down all the 12 identities 
fl4.16p - fl4.18p in the compressed form 


(-i)9+’9?(«)9j+i(i.)-9yi(«)9yi(«) 

9„+i(2t.)9yi(0) 

9„+i(2t.)9yi(0) 


9j+i(2u)9„+,(0)9y,(0), 

(4.19a) 

-9,+i(2«)9,+i(0)97i(0). 

(4.19b) 

c,(«) + (-ireiiv . 

(4.19c) 

(-i)’+'9y,(«) + (-i)'’+‘97,(«). 

(4.19d) 


where in fl4.19ap . fl4.19bp . and fl4.19dp the indices a, /5 ,7 are assumed to be any cyclic permu¬ 
tation of {1,2,3} and in fl4.19cl) a = 1,2,3. Together with fl4.15p . the system fl4.19p yields 
the complete set of duplication formulae. 
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